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Abstract 

We present a sample of edge-on spiral galaxies both of early and late types. 
The sample consists of 175 galaxies in the ET s -filter, 169 galaxies in the i^-fllter 
and 165 galaxies in the J-filter. Bulge and disc decompositions of each galaxy 
image, taken from the Two Micron All Sky Survey (2MASS), were performed. We 
discuss several scaling relations for bulges and discs which indicate a tight link 
between their formation and evolution. We show that galaxies with bulges fitted 
by the Sersic index n < 2 (pseudobulges) have quite different distributions of their 
structural parameters than galaxies with n ^ 2 bulges (classical bulges). First 
of all, the distribution of the apparent bulge axis ratio % f° r the subsample with 
n < 2 can be attributed to triaxial, nearly prolate bulges, while n ^ 2 bulges seem 
to be oblate spheroids with moderate flattening. Secondly, the Photometric Plane 
of the sample bulges is not flat and has a prominent curvature towards small values 
of n. Thirdly, despite of the existence of a clear relation between the flattening of 
stellar discs h/z$ and the relative mass of a spherical component, the distributions 
over both parameters are quite different for galaxies possesing classical bulges and 
pseudobulges. 

Keywords: Galaxies: infrared — photometry — spiral edge-on galaxies — fundamen- 
tal parameters. 

1 Introduction 

In recent decades a number of studies concerning structural properties of galaxies has 
been performed. The results of decompositions of near face-on disc dominated galaxies as 
well as of early-type galaxies yeilded a lot of correlations between structural parameters 
of discs and bulges [U EJ EJ S]. 

Among all galaxies edge-on ones are of great interest because they provide a unique 
possibility to obtain information about the vertical structure of a disc and a bulge. The 
disc scaleheight z Q together with the disc scalelength h determine the relative thickness 
of a stellar disc for each galaxy. This ratio as well as the bulge flattening % (the ratio of 
isophote semi-axes) may put constraints on the disc and bulge formation processes and 
their secular evolution. 

There are some difficulties in studying edge-on galaxies and obtaining their structural 
parameters because of the dust extinction. The possible solution is to analyse observa- 
tions in the NIR passbands where the extinction is drastically reduced. A detailed study 
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Table 1: Morphological types of the sample galaxies and their fractions 



Type 


Number 


Percentages 


Type 


Number 


Percentages 


SO, SO-a 


41 


23.4 


Sbc 


21 


12.0 


Sa 


11 


6.3 


Sc 


36 


20.5 


Sab 


15 


8.6 


Scd 


4 


2.3 


Sb 


47 


26.9 


Total 


175 


100 



of the discs of edge-on galaxies is described in many papers [5j EJ [TJ El El QUI [TTJ. Almost 
all these papers dealt with bulgeless or late types edge-on galaxies and did not focus on 
the properties of bulges. Our study is intended to investigate structural parameters of 
edge-on galaxies both of early and late types and to join analysis of their bulges and discs 
in the near-infrared bands. 

For our purpose we used the 2MASS-selected Flat Galaxy Catalog (2MFGC) rj2] as 
a source of objects and the 2MASS survey [32] as a source of IR image data. We also 
used the Revised Flat Galaxies Catalog (RFGC) [H] to add late-type galaxies to our 
sample. Our celection criteria are presented in [15]. The absolute numbers of galaxies of 
various morphological types (taken from the LEDA) and their percentages are listed in 
the Table [U 

The distances of the galaxies are taken from NASA/IPAC Extragalactic Database 
(NED) with the Hubble constant H = 73 km s _1 /Mpc, ^matter = 0.27, and fi vacuum = 
0.73). To investigate some kinematical and dynamical properties of the sample galaxies 
we added information about maximum rotation velocity for each galaxy if it was presented 
in the LEDA data base. The sample consists of 175 galaxies in the i^ s -filter, 169 galaxies 
in the if-filter and 165 galaxies in the J-filter. The full list of the sample galaxies with 
the basic information about them can be found in [15]. 

According to the V/V max test [TF] our full sample is incomplete but the subsample of 
galaxies with angular radius r ^ 60" is complete (V/V max = 0.49 ±0.03). This subsample 
consists of 92 galaxies (47 early-type galaxies and 45 late-type galaxies). 

2 Two-dimensional bulge/disc decomposition 

We used BUDDA (Bulge/Disc Decomposition Analysis) v2.1 code [18] to perform two- 
dimensional decomposition of galactic images taken from the database of the All-Sky 
Release Survey Atlas. 

The photometric model consists of two major components: a bulge and a disc. The 
disc is assumed to be axisymmetric, transparent and is represented by an exponential 
distribution of the luminosity density. It can be described by the face-on central surface 
brightness Sq,<i (in mag per arcsec 2 ), the scalelength h and the 'isothermal' scaleheight 
z . 

For the bulge surface brightness profile the Sersic law is adopted [16J. The main 
parameters of this profile are: the central effective surface brightness /io,b expressed in 
mag per arcsec 2 , the effective radius of the bulge r e ^ and the Sersic index n, defining the 
shape of the profile. We also assume that bulge isophotes are described by ellipses with 
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the flatness % — b/a, where b is the semi- minor axis and a is the semi- major axis of an 
isophote. The Sersic index of n = 4 represents the de Vaucouleurs profile which is very 
popular to describe the surface brightness distribution of bright elliptical galaxies and of 
bulges in early-type spirals. The Sersic index of n = 1 represents the exponential profile 
of bulges in late-type spirals, of galactic discs and of dwarf elliptical galaxies (see e.g. 

na m m m)- 

To check the reliability of our 2D decompositions, we compared our results with the 
results obtained by other authors. In the K s filter, our sample has 30 galaxies in common 
with the sample in [7] and 14 galaxies that have been investigated in j5]. All mentioned 
authors used different decomposition procedures but the results of comparison are in 
good agreement [T5] . 

3 New evidences for different origin of bulges and 
pseudobulges 

Statistical correlations between structural parameters give us a clue for understanding 
the formation and evolution of spiral galaxies. Here we present several correlations of the 
main parameters of the galaxies and some new evidences in favour of existence of two quite 
different populations of bulges in spiral galaxies — classical bulges and pseudobulges. 

For edge-on galaxies morphological types are very subjective because the spiral arms 
are not seen in this case. That is why we distinguish the galaxies in our sample using the 
Sersic index n because there is a good correlation between n and the ratio of bulge and 
disc luminosities B/D that is known to be a key parameter for classifying spiral galaxies 
on the Hubble sequence JTJ [3] 

The distributions over the Sersic indices of galaxies in three passbands are shown 
in Fig. [T] (the left plot). The distributions demonstrate a weak bimodality which may 
reflect the existence of two families of bulges: bulges with n ^ 2 and bulges with n < 2. 
This bimodality is known for some galaxy samples (e.g. [22]) and probably is real for 
our sample. As was reported earlier (see e.g. [221 [23]) such a bimodality correlates with 
morphological type of the bulge. The classical bulges have n ^ 2 and the so-called 
pseudobulges have n < 2. The classical bulges appear to be similar to elliptical galaxies 
and have similar properties. It is assumed that these systems were built via minor and 
major merging. The galaxies with the pseudobulges are thought to be formed via disc 
instabilities and secular evolution and have disc-like apparent flattening. 



3.1 The intrincic shape of bulges 

In Fig. [T] (the right plot) we demonstrate the distributions of model bulge axis ratio 
for edge-on view. In the case of edge-on galaxies this parameter describes directly the 
bulge flatness and the intrinsic 3D structure of bulges if they are assumed to be oblate 
spheroids. The median value of this parameter is ~ 0.63, independently of the band. This 
is in good agreement with [SUES]. Thus, in general bulges are definitely nonspherical and 
flattened. However the question about 3D shapes of bulges and their possible triaxiality 
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Figure 1: Distributions of the sample galaxies over the Sersic index of bulges (left) and the 
model bulge axis ratio in the iT^-band (right). At the right panel the thin line indicates the 
distribution for the whole sample, the solid line corresponds to the subsample of bulges with 
n 2 and the dashed line shows the distribution for the subsample of bulges with n < 2. 



are still open, despite of the fact that its solution is thought to be crucial for testing 
different scenarios of galaxy formation. 

Our sample of edge-on galaxies allows us to distinguish clearly the difference in the 
bulge equatorial ellipticities for galaxies with classical bulges and pseudobulges. As one 
can see in Fig. [1] (the right plot) for classical bulges the distribution of has a rather 
narrow peak at % ~ 0.65. This may reflect the fact that the bulges are nearly oblate 
spheroids with moderate flattening. The distribution over % for pseudobulges is very 
wide, spreading from very flat bulges up to nearly spherical ones. Such a distribution 
may be attributed to definitely triaxial, near prolate bulges that are seen from different 
projections — along the major axis and perpendicular to it. But we can not exclude that 
triaxial shape of bulges in late-type galaxies may hide the presence of bars that went 
thick in the vertical direction during secular evolution. 

Our findings conradict the results obtained in |26j. These authors recovered the PDF 
for the sample of 148 unbarred SO-Sb galaxies at intermediate inclination angles and 
concluded that bulges with n < 2 show a larger fraction of oblate axisymmetric bulges, 
a smaller fraction of triaxial bulges, and fewer prolate axisymmetric bulges with respect 
to bulges with n ^ 2. 

3.2 Photometric Plane and Sersic index 

The difference between bulges and pseudobulges arises while building the Photometric 
Plane (PhP) for them. This plane represents a tight correlation of the Sersic index n 
with the central surface brightness yU ,b an d the effective radius of a bulge r e b J27J 12H] • 
Performing the least-squared fit of an expession lg n = a lg r e + b /i b + c we obtained a 
surprising result. It occurred that the galaxies with Ign ^ 0.2 (n ^ 1.58) lying on the 
PhP, that is built for these galaxies, has a small scatter, but the galaxies with Ign ^ 0.1 
(n ^ 1.26) do not correspond to this plane and form their own plane (Fig. [2]). The 
main difference between our sample and the samples of the above-mentioned works is the 
range of the shape parameters n. Our sample contains a substantial amount of bulges 
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Figure 2: The Photometric Plane built for bulges with lgn ^ 0.2 of the complete subsample 
in the i^-band. The points corresponding to bulges with lgn ^ 0.1 do not lie in this plane. 
The model merger remnants of disc galaxies are shown with open circles and triangles [29J. The 
model data are shifted by a constant along the horizontal axis. 

with lgn ^ 0.1. For the sample in [27J the low boundary for n lies at lgn ^ 0.2, that is 
why the curvature of the Photometric Plane towards small values of n for bulges was not 
noticed earlier. 

It is not clear whether two distinct parts of the curve in Fig. [2] reflect the different 
origin of bulges lying in the corresponding regions. However, the collisionless merger 
remnants of disc galaxies obtained in iV-body simulations and fitted by a Sersic profile 
fairly well reproduce the slope of the PhP in the region of lgn ^ 0.2 [29]. We put on 
the data from [29] in our Fig. [2] and found the numerical data to be consistent with our 
observational ones. The curvature of the PhP towards small values of n may reflect the 
quite different nature of such bulges, formed, for example, via secular evolution of discs. 

The found curvature is not the feature of our decomposition and the presence of a large 
number of galaxies with small angular radii which may give a confusing decomposition. 
We constructed the PhP only for galaxies from the complete sample, possesing systems 
with r ^ 60". The curvature was not noticed earlier because of a lack of substantial 
number of late-type galaxies in the samples that were studied. We constructed the 
Photometric Plane based on the samples presented in [301 [31] • The first sample contains 
121 face-on and moderately inclined late-type spirals with derived structural parameters 
of bulges and discs in the if-band. Fig. [3] (the left panel) shows the Photometric Plane 
for the bulges from [30] together with our data. The plot displays the same feature as 
Fig. [2j a very tight correlation for classical bulges and a fairy large scatter of points 
for pseudobulges. Surprisingly, both samples follow the same relation in the region of 
small values of n. The sample, presented in |3T], consists of around 1000 galaxies with 
determined photometric parameters in the i-band. The most close passband in our sample 
is the J-band. We display the best fitting relation for our data in Fig. [3] (right) and added 
the point from the sample [31J. There is a shift between two relations which is due to 
the difference in the passbands. But both samples show a prominent curvature towards 
small values of n. 
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Figure 3: The Photometric Plane built for bulges with lgra ^ 0.2. Left — the data from 
|30j (open squares) in comparison with our data (filled circles). Both data are related to the 
f/-band. Right — the data from [5T] (open squares, the i-band) in comparison with our data 
(filled circles, the J-band). 

3.3 The dark halo and disc flattening 

The disc flattening can be measured directly only for edge-on galaxies. It is commonly 
expressed as the ratio of disc scaleheight z to disc scalelength h. This ratio shows a 
weak trend with Hubble type [5j. Zasov et al. [10] were the first to found that the 
disc flattening shows a definite trend with the ratio of dynamical mass M tot to the disc 
luminosity provided the dynamical mass is defined as the total mass enclosed within 
the sphere of radius that is equal to four disc scalelengths. These authors suggested 
that discs are marginally stable against the growth of perturbations in their planes and 
bending perturbations and implied the linear relation between the ratio h/zo and the 
ratio of dynamical mass (or dark halo mass) to disc luminosity. It is valid at least for 
bulgeless galaxies. 

The found correlation means that relatively thinner discs in bulgeless galaxies tend 
to be embedded in more massive dark haloes. Zasov et al. [10] verified the linear trend, 
obtained analytically, for two different samples of edge-on bulgeless galaxies with known 
structural parameters in the R and K s bands, using the HI line width W^q as a measure 
of dynamical mass. For our sample, which contains a lot of galaxies with massive bulges, 
we performed the same analysis. 

We chose the ratio M to t/Md to demonstrate a trend for the disc flattening with the 
relative mass of a spherical component (Fig. H]). We adopted mass-to-light ratios fj = 1.5, 
fn = 1-0, and fx s = 0.8 M Q /L Q [22] to get an estimate of the disc mass M d . The 
scatter of points in Fig. @] is rather large in contrast to the results displayed in [ID] . 
It may be caused by uncertainties in velocities and mass-to-light ratio determinations, 
decomposition errors and disc relaxation processes, that can thicken a disc just above the 
threshold for marginal stability. But some of this scatter must be real. Sotnikova and 
Rodionov [33] concluded that the presence of a compact bulge is enough to suppress the 
bending instability that makes the disc thickness increase. A series of iV-body simulations 
with the same total mass of a spherical component (dark halo + bulge) were performed. 
The final disc thickness was found to be much smaller in the simulations, where a dense 
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Figure 4: The ratio of h/zo as a function of the ratio of dynamical mass to disc mass. The 
open squares correspond to the subsample of bulges with n ^ 2, the filled squares refer to the 
subsample of bulges with n < 2. Specially marked points, which deviate from the main relation, 
refer to the galaxies with peculiarities in their structure (see |15|). 

bulge is present, than in the simulations with bulgeless systems. The results of iV-body 
simulations of discs starting from an unstable state were summarized in [31]. The ratio 
z /h versus (M h + M b )/M d (M h is the halo mass within the sphere of radius Ah and M b is 
the bulge mass) was plotted and it was shown that there is a clear scatter in this relation, 
in spite of the same model mass for a spherical component (Mh + Mb). 

Fig. H] shows no difference between galaxies with bulges and pseudobulges but there 
is a hint of bimodality in the distribution of our sample galaxies over the ratio of the 
dynamical mass to the stellar mass M* = Mj + Mb (Fig. [5j the left plot) and over the 
ratio h/z (Fig. [5j the right plot). This may reflect the presence of two different families 
of galaxies with different bulges in our sample. Galaxies with pseudobulges posess the 
large fraction of a dark component and their discs are relatively thiner. 



4 Summary and conclusions 

We constructed the large sample of edge-on spiral galaxies with performed 2D bulge/disc 
decomposition of 2MASS galaxies images in J-, H- and A's-passband using the BUDDA 
v2.1 package. We extracted global bulge and disc parameters for all objects. 

Analysing correlations between structural parameters, including those that describe 
the vertical structure, we found new evidences in favour of existence of two quite different 
populations of bulges in spiral galaxies — classical bulges and pseudobulges. Our results 
can be summarized as follows. 

1. The bulge flattening % divides the sample into two different families — triaxial, 
nearly prolate bulges and close to oblate bulges with moderate flattening. The 
Sersic index threshold n ~ 2 can be used to identify these two bulge types. 
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Figure 5: Distributions of the sample galaxies over the ratio of dynamical mass to stellar mass 
of a galaxy (left) and the ratio of h/zo (right). The dotted line indicates the distribution for 
the whole sample, the dashed line corresponds to the subsample of bulges with n #s 2 and the 
solid line shows the distribution for the subsample of bulges with n < 2. 

2. The most surprising result arises from the investigation of the Photometric Plane of 
sample bulges. The bulges with n ^ 2 populate a narrow strip in their Photometric 
Plane. However, there is a difference in behavior of this plane for bulges with n ^ 2 
and n < 2. The plane is not flat and has a prominent curvature towards small 
values of n. This result may be due to the physical distinction between classical 
bulges and pseudobulges. 

3. There is a correlation between the relative thickness of stellar discs h/z and the 
relative mass of a spherical component, including a dark halo. This correlation 
was known previously for bulgeiess galaxies [10] or samples with predominantly 
bulgeiess galaxies [9] and was argued to arise from marginally stable discs. Our 
sample is much larger than the samples in 0, [TU] and more reliable for statistical 
analysis. Here, we confirm the correlation under discussion. What is more, we did 
this not for bulgeiess galaxies but for galaxies with massive bulges and concluded 
that the total mass contained in spherical components may be one of the factors 
that determines the final steady state disc thickness. 
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